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1. Introduction 

Wigner and Husimi functions are important tools both in quantum mechanics and signal 
analysis. They both provide a phase-space description of wave functions and signals. 

Wigner provides a formulation of quantum mechnaics that seems , at first sight, very 
close to transport equations for dynamics of assembly of classical particles. However an 
important and inconvenient feature of Wigner functions is the fact that they are not 
pointwise nonnegative. This important difference with classical densities causes not only 
a problem of interpretation, it also creates technical problems. 

The Husimi function is usually seen as the result of an appropriate smoothing the 
Wigner function, and this smoothing is enough to make it positive. The simplest way of 
seeing this is to observe that the Husimi function is nothing but the square modulus of the 
scalar product of the original function with a "gaboret" , namely a (Gaussian) coherent 
state. 

In the case where the original function satisfies a partial differential equation, e.g. a 
Schrodinger equation, Wigner and usual Schrdodinger representation s have been known 
to be, modulo a global phase, equivalent. Given a Schrodinger equation, there is a corre- 
sponding Wigner equation (see e.g. [7]) which is in general a pseudodifferential equation 
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and whose limit, as the Planck constant h vanishes, is the classical Liouville equation. Al- 
though Wigner and Liouville equations are "close" for small values of h, the non-positivity 
of the Wigner function forbids to interpret it as a classical density, and therefore the 
Wigner equation as a perturbation of the Liouville one inside the paradigm of classical 
mechanics. It seems therefore natural to somehow try to overcome this difficulty by de- 
scribing the quantum system in the Husimi picture. However, by doing that, we encounter 
immediately the difficulty that the benefit of positivity is balanced by a dynamics less easy 
to derive. This is due to the fact that smoothing is hard to "effectively" invert. 

Recently, [1] , the equation that is satisfied by the Husimi function of a function, solution 
of a Schrodinger equation, has been derived. This "Husimi equation" involves in general 
analytic continuation in the complex domain, and can be expressed without complexifi- 
cation of the arguments in case of analytic Hamiltonians, for example when the original 
Schrodinger operator is a differential operator with polynomial coefficients. 

In this paper we will derive the corresponding equation for the affinc Husimi function. 
After a review of the standard Weyl- Wigner situation and the main results of [T] in Section 
[2j we will define define "Weyl" quantization in the "ax+b" group setting in Section [H and 
show how, using continuous wavelet (affine coherent states) one can define positive Husimi 
functions in Section 01 we will then derive the corresponding Husimi equation and express 
our main theorem Theorem l5.2l in Section [5] and suggest possible generalizations in section 

m 

2. The usual Weyl case 

In this section we briefly review the standard smoothed Wigner results (see [1] for 
details). 

The fundamental equation of quantum mechanics, namely the Schrodinger equation, is a 
partial differential equation involving a self-adjoint operator on a Hilbert space, generating 
a unitary flow. Usually the Hilbert space is L 2 (W l ) and the unknown is the so-called wave- 
function. This fact makes quantum mechnaics different from classical one by two facts : 
the fundamental equation is a PDE and there is no, a priori, reference to a phase-space, 
despite this it has been recognized since the eraly days of quantum thery that a phase-space 
description of teh wave-function is possible and useful. This fact has to be put together 
by the time-frequency methods in signal analysis. One of the best ways to illustrate this 
common feature is to introduce the so-called Wigner- Ville transform. 

On M. n the Wigner function associated to an L 2 (R") function t/j is by definition: 



(2.1) W i ,{x,£)= e- 2i ^ x <P(x + h$)1>(x - h%)dy 
and its smoothed version is given by 

(2.2) W^(x,0 = (V2/fla x a i ) n e ~^~ h ^W^V, ?)<&(!?. 

It is well know that if ip satisfies a Schrodinger type equation of the type: 

(2.3) ihd t tfj — h(x 7 ~ihV x )ip, 
the Wigner transform satisfies the following equation: 

(2.4) OtW^ = 23? (ih(x ~ ~ hV £ , £ + l -hV x )Wip 
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In what precedes h(x, —ihSJ x ) is a pseudodifferential operator given by the Weyl calculus: 
(2.5) h(x, -ihV x )u{x) = h~ n { h((x + y)/2, ^)e' ix - y)/h u(y)dy 



and h(x — |^V^,^ + |W X ) is pseudodifferential operator on L 2 (R 2n ) obtained the same 
way. 

In 1 was proven the following theorem: 

Theorem 2.1. Let ip satisfies fH3J) with (for simplicity) h € <S(M 2 ") (the Schwartz class). 
Then W^t satisfies the following: 

(2.6) d t W r = 2K (JCW^A 

where 

Cw(x,0 := J h(S,T)e^ s * +T V-^ s2 +°t T2 )w ( x + h T + i ^ S ^-n^^\ dSalT. 

(the possibility of extending the argument of w is clear from the fact that the smoothed 
Wigner function is analytic on the whole complex plane). Here h is the Fourier transform 
ofh. 



There are several ways of "formally" deriving ()2.6[) . The first one consists in smoothing 
and unsmoothing the equation (J2.4I) and apply Egorov theorem associated to heat equation 
(i.e. free evolution at complex time). 

The second one consists in deriving the formula from the smoothed product itself. 

Finally, in the case a x x Ok — 1 (and only in this case) we can derive (|2.6|) by using the 
decomposition of identity on coherent states formula. This is the method we will use for 
the wavelet case. 

Let us remark mention that, in the case where the symbol h is analytic (more precisely 
is the Fourier transform of a compact supported function) one can write th evolution 
equation for W^ as: 



d t W r =2^(c a W r ) 



where C a has Weyl symbol: 



/ K - icr 2 X X + ialK 
L(x,t:X,K):=h(x-h -^,£-ft -*— 



3. The affine Weyl quantization 

Inspired by [5] affine Weyl quantization has been introduced in ([10). We set the 
following . Let on L 2 (R + , dx) , a > 

U(a,b)(p(x) :— a l ' 2 e~~ ip{ax) 

U(a, b)U(a', b') = U(aa', ab' + b). 
Let Itp(x) := (p(—x) In analogy with the Weyl case we define: 

V{a,b) :=U(a,b)IU(a,b)-\ 
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and the Weyl quantization of W(a, b) the operator: 

(3.1) W= f W(a,b)V(a,b) dC "' h 



2a 2 h 



With the condition that W £ L 1 (M+ x R, ^) (|3TT|) defines clearly a (non uniformly in 
K) bounded operator on L 2 (M. + ). 

Remark 3.1. We will not enter here in boundness condition uniform in h (Calderon- 
Vaillancourt type results). We refer to [2] for an extensive study. 

Lemma 3.2. ifW is given by H3.1\) , then its integral kernel is (note the analogy with the 
standard Weyl quantization): 

(3.2) W ( xy ) = w(-=,y^z), 

\/xy 



W{a,c) = -^= I W{a,b)e~ l ^db. 
\j2irh . 



where W(a,c) is the h-Fourier transform in the second variable: 

(a,c) 

In particular J R+xK V(a, b)^^ = Identity. 
Proof. We first notice that U(a, by 1 = U(±,-%). therefore: 

V(a, b)xj;{x) = ^(-Ly^-^, 
a z x 

from which we deduce that the integral kernel of V(a, b) is 5{y ^-^ e l (v- x ) b / h from which 

we get ([3T2]) . D 

An easy computation shows the following 

Lemma 3.3. Let C(a) the operator defined on L 2 (M. + ) by: 

C(a)ip(x) ■- - { ax + —J fix)- 

We have: 

(3.3) W(a, b) = Tr[WV(a, b)C(a)} . 

Remark 3.4. the extra factor C(a), not present in the standard Weyl quantization, is due 
to the non-unimodularity of the u ax + &" group. 

We define the affine-Wigner function of a function ip £ L 2 (R + ) as the Weyl symbol of 
the orthogonal projector on ip, divided by the Planck constant. 

Definition 3.5. 

W^(a,b) := -Tr[WV{a,b)C(a)} = (ip,V(a,b)C(a)il>}. 



We get immediately the following expression: 



W^,(a,b) = — / (aaH ) i/;(x)e 

2hJ R + \ ax J 



'■"""fcVijfcl*.. 



and 
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A study per se of the semiclassical limit of W x p(a,b) has not, at our knowledge, done 
in the spirit of weak convergence results for the usual case [7J |5] and will be done in [2] . 
In this paper we will study "directly" its smoothed version since, as in the standard Weyl 
case there is no reason for which the affine-Wigner function should be positive. 

4. Affine Husimi 

Coherents states have a long history in quantum mechanics (they were introduced by 
Schodinegr in 1926), and also have got a more recent interest in signal analysis through the 
so-called Gaborets. The afiine setting ihas produced the wavelets, and in qm corerspond 
to a situation where the undelying phase-psace if the upeer half plane. 

We define now the set of coherent states (continuous wavelet): 

Definition 4.1. 

fa,b{ x ) = C{h)a 1 ' h+ ^x l / h e-^ a -^ x / h = U(a,b)(fi,o(x). 



where C(K) = \ jb-z ", ,■, is a normalization coefficient insuring that 

1 1 </?a,b || L 2 (It+)- 

We first remark that ip a ^ decreases exponentially at +oo and as an increasing power 
as h —> at the origin. Ir results from general "wavelet" theory (and a straightforward 
computation) the following decomposition of identity: 

i a \ f i dadb 

(4.1) / \y a ,b >< <p a ,b\ — z- = Identity 



1) 






/ \(fia,b >< <fia,b\ 

Jlxl+ 


dadb 
ah 


Wc 


i define 


now 


the wavelet-Husimi function: 





Definition 4.2. The Husimi function of an operator W of integral kernel W(x, y) on 
L(R+) is s 

W(a,b) •■= {iPa ' b ' WiPa ' b) ■ 
h 

where < ., . > is the scalar product on L 2 (M + ). 

Let us first remark that, thanks to (|4.ip we get immediately: 

f — dadb 

(4.2) / W = TrW. 

Jrxm+ a 

In the case where W is the orthogonal projector W^, on a function if) £ L 2 (R + ) (pure 
state) we get that 



(4.3) 


W^(a,b) := C(h) 2 \ / x 1/n e' ia+ib)x/h ip(x)dx 


and: 




(4.4) 


Jlxl+ a 



Let us remark finally that, by construction, W^ > and W > if W > as an operator. 
It is easy to compute the wavelet-Husimi function out of the affine-Wigner one defined 
earlier. 
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Proposition 4.3. Let us denote wm{s,o) the Mellin transform in a of w(a,b) (see the 

Appendix) . We have: 

(4.5) 

_ ol+2/fi r / 2 , /L t\2\ f-TT- 1 

^(M)= tr , 2x . / \ ) T( S /2-l/h-lfw M (s,0dsdt; 

hl \h + L > J(l/2+iR)xR V " / 

Therefore ifW>0 as an operator (for example if W — \ip >< ip\ for some vector ip) and 
if W(a, b) has a weak-limit the resulting measure is positive. 

Proof. We first write that 

((Pab,Wipab) 



W(a,b) 



h 

Therefore, since the integral kernel of W is W(—l=,x — y) we get, using the inverse of the 
Mellin transform given in the Appendix, 

p pl/2-\-ioo 

W(a,b) = C{h) 2 / / {xy) s/2+l/h e- {a+l{b - i))x - {a - l{b -^ )v w M {s^)dxdyd^ds 

J(R+) 2 Jl/2-ioo 

using the equality: 

A- fc r(fc) = / u k - l er uy du 

Jr+ 

we get the (j4.5[) . It is easy to check that, as h — > 0, the r.h.s. of (|4.5p tends to w(a, b) if 
this latter has a limit wq as a measure. Therefore, since W > we get that wq must be 
positive. 

□ 

5. The result 

We can now state the main result of this paper. 

Theorem 5.1. Let H be an operator of affine-Weyl symbol h in L 1 (R + x M, ^pr). Let 
%l> be the solution of the Schrodinger equation; 

(5.1) ihdti! t =H^ t . 

Then the wavelet-Husimi transform of ip* , Wl satisfies the following equation: 

(5.2) d t W r (a, b) = ^ (J $(a, a; b, (3)W r (a + ^^,b + ^±^)dad[3 

where 

2 1 + 2 / fi f f (a + ijP - Q(a + 2a- »(6 + 2/3 - fl \ * ~^~ 1 

$(a,a;6,/3) = 2 / z x 

ftl(ft + 1) 7(1/2+-*) xr V n J 

a + a + i(b + /3) 2 ... 

T(s 2-1 h-1) w M {s,Qdsd£ 

a — ib 

and 9z is the imaginary part of z. 

The same property holds for a non-pure state W satisfying ihdtW = [H, W ] by re- 
placing ip* by W* . 
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Theorem [lO] is actually a direct consequence of the Theorem l5.2l below. Indeed we note 
first that the equation (|5.1[) can be rewritten as: 

a t \i> t ><i> t \ = ^[H,\ii> t ><ii>% 

in 
Therefore 

dt(tp(a,b),W X^WiaM)) = (V(a,6),^[-H'>IV'*><V' i |]^( a ,6)> 

= -TT «¥V,6)#|V>* >< ^*IV(o,6)> ~ (VW)!^* >< V»*l<HV{o,6)» 
= ~$ (\V (a , 6 )F|^ X V*k(a,6)» 

Theorem 5.2. 

' — ' /" Q; — %3 3 ~\~ tot 

(5.3) W^(a, b)= *(o, a; 6, ^(a + — ^-A 6 + ^— )dad/3, 

w/iere 

t („, 0iM) . *£L/ A. + ^-0(» + 2a-,( t + 2^-m *-*-' x 

fii VI + 1 i J(l/2+iR)xK V ft / 

a + a + i(b + (3) 2 , - , , t 

a — lb 
where WM(s,b) is the Mellin transform in the first variable. 

Proof, the proof will use as one of the main ingredients the following Lemma. 

Lemma 5.3. The wavelet-Husimi transform W^(a,b) has an analytic continuation which 
is given by the following formula; 



(5.4) Wj,(a + a + iP,b + P-ia) = J — — (<P( a +2a,b+2P)i W( Pa,b)) 

Proof D 

We now use the decomposition of the identity by afhne coherent states (wavelets) that 
is: 

._ _. f . dadb 

(5.5) / \(p( a> b) >< <P( a ,b)\ — z- = Identity. 

plugging (J5.5I) in: 

dad/3 

we get the result thanks to the following Lemma, easy extension of Proposition ^. 31 
Lemma 5.4. 

{<P(a'P),W<p [a , b) ) = 2 / ~ X 

M Vft + J J «/(l/2+*R)xK V ft / 

T(s/2-l/h- l) 2 w M (s,£)dsd£ 

n 



(5.6) (p( 0>6 ),.ffWV(a,&)) = / (^( a ,&),-H r WV(a+2a,&+2 J 8)) (<£(a+2 Q , 6+2/3), W^fj,!,)) 
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6. Possible generalizations 



One of the interest of coherent states is the fact that, since they are purely local, they can 
be defined in situation where there is no corresponding Weyl calculus, in particular in case 
of equations on manifolds [9] . They even are the building blocks of semicalssical methods 
in the case where the phase-space is not even of the form T*A4 (cotangent bundle), and 
is for example, a Kahlerain manifold. 

We believe that it would be possible to derive the Husimi calculus in theses situations, 
and that it should have useful applications to PDEs theory (see [3j [4] for applications of 
afiinc Wigner functions in signal analysis). 

Appendix A. The Mellin transform 

In this short appendix we recall some very well known facts about the Mellin transform. 
To a function w defined on the positive axis we associate its Mellin transform wm 
defined through: 



M') 



/ x s w9x) — . 



The transform w(x) — » (2n) 1 / 2 w(l/2 + it) is a unitary transform between L 2 (R + ) and 
L 2 (M) and, therefore, the Mellin transform can be inverted thanks to the formula: 

-i />l/2+ioo 

w(x) = -— / x~ s w M (s)ds. 

2m Jl/2-ioo 

Let us note the link between Mellin transform and Laplace L and Fourier ^transforms: 

Fw(0 = (wo (-log)) M (i0, 
Lw(t)) = (too (-\og)) M (i]). 
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